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We show that a weakly absorbing, strongly scattering (white) medium can be made very strongly 
absorbing at any frequency within its strong-scattering bandwidth by optimizing the input electro- 
magnetic field. For uniform absorption, results from random matrix theory imply that the reflectiv- 
ity of the medium can be suppressed by a factor ~ (ia/i)N~'^ , where N is the number of incident 
channels and i,ia are the elastic and absorption mean free paths respectively. It is thus possible 
to increase absorption from a few percent to > 99%. For a localized weak absorber buried in a 
non- absorbing scattering medium, we find a large but bounded enhancement. 

PACS numbers: 42.25.Bs, 42.25.-p, 05.60.Cd 



The absorption properties of a medium at a given fre- 
quency depend not only on the atomic and molecular 
constituents of the medium, but also on the coherence 
properties of the incident radiation. This was empha- 
sized by recent work demonstrating the phenomenon of 
coherent perfect absorption (CPA), in which a cavity con- 
taining a weakly absorbing medium completely absorbs 
an appropriately-chosen input field; this is achieved by 
choosing the field to be the time-reverse of the lasing 
mode which the cavity would emit, if the loss medium 
were replaced by a gain medium of equal strength [U |2] . 
Other input fields, differing only in the phase relation- 
ships of the incident waves, are only weakly absorbed, 
as are incoherent inputs. Like a laser, CPA only occurs 
at discrete frequencies in a given system. This raises a 
question: under what conditions, if any, is it possible to 
control absorption over a continuous and large frequency 
range by optimizing the input fields at each frequency? 

A related question, which has recently been extensively 
studied, concerns the conditions under which light can 
be transmitted through an opaque, lossless^ multiple- 
scattering medium. In the diffusive regime, the trans- 
mission probability decays as where £ is the elastic 
transport mean free path and L is the sample length; 
hence a sample with L / is highly refiecting at all fre- 
quencies. Nonetheless, it has been known for some time 
that at any frequency there do exist coherent superposi- 
tions of input fields which penetrate much further than 
the mean free path, and hence may transmit with high 
probability through the sample [3-6 . This effect has 
been demonstrated experimentally in the past few years, 
as practical methods were found for determining the op- 
timal input fields without detailed a priori knowledge of 
the scattering configuration [71 |8l . This suggests that if 
a multiple-scattering medium contains some material ab- 
sorption, then even if that absorption is weak (i.e., the 
medium normally appears "white" ) , it should be possible 
to substantially enhance the effective absorption at any 
frequency within a wide and continuous frequency range, 
by tuning the input fields to penetrate deeply enough 
that the photon path length exceeds the large but finite 
absorption length. In this Letter, we provide a theoreti- 
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FIG. 1: (color online) (a) Model of an absorptive multiple- 
scattering system. A waveguide with dielectric constant is 
connected to a region of length L with dielectric nl -\- Se-\- ir/, 
where 6e{x) is a real white-noise disorder term drawn from 
a uniform distribution over [—do, do], and 77 is the material 
absorption, (b) Variation of reflection eigenvalues {Rn} with 
frequency k for a fixed disorder realization, with no = ni = 
1.5, N = 80 scattering channels, mean free path i = 0.05L, 
and uniform 77 with ballistic absorption length ia = 13.3L. 
Dashed lines show the ensemble averages of the mean and 
smallest reflectivity, R and Ri . (c) Semi- log plot of Ri , which 
fluctuates, but is ^ 10 ~^ over a continuous frequency range. 

cal and numerical demonstration of this effect, which we 
refer to as coherently enhanced absorption (CEA). It is in 
general distinct from CPA, which follows rigorously from 
time-reversal symmetry and is not specific to multiple 
scattering media (but requires fine-tuning of the sample 
parameters and operating frequency). 

Linear scattering from a medium is described by the 
scattering matrix (S-matrix), which relates incident and 
scattered amplitudes in the basis of asymptotic free solu- 
tions of the wave equation. The S-matrix depends on the 
scattering geometry and material properties, as well as 
the frequency of the input fields, but does not depend on 
the fields themselves. Given a vector |?/^in) whose com- 
ponents are the intensity-normalized input amplitudes 
in each of the N scattering channels, the output ampli- 
tudes are IV^out) = and the total output intensity 
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is (V^outlV^out) = {i^in\S^ S\ilJin)' For a lossless system, this 
quantity is always unity since 5* is unitary; in the pres- 
ence of absorption, it is smaller than one [9 . 

We will focus on the case of scattering from one surface 
of an optically dense and weakly absorbing medium, in 
which all incident fields are either reflected or absorbed, 
so that the S-matrix coincides with the reflection matrix 
r. Denote the eigenvalues of the Hermitian matrix 
by Rn^ where < Ri < - ■ < Rn < 1- The sample- 
specific "mean reflectivity" R = (l/N) Rn is the nor- 
malized output power when the inputs in each channel 
are independent (e.g. having no coherent phase relation- 
ship); the moments of this quantity were studied in ear- 
lier works on reflection from random absorbing media 
[lOl[TT]. In studying CEA, the relevant physical quantity 
is the smallest reflection eigenvalue, Ri. By the vari- 
ational principle, (V^in^^^l'^in) ^ for any intensity- 
normalized IV^in)- Hence, 1 — Ri is the maximum amount 
of absorption that can be achieved, and its eigenvector 
gives the maximally- absorbed set of input amplitudes. In 
an experiment, if noise sources rule out measuring reflec- 
tivities below a value R' > R^^ any input waveform in 
the eigenspace spanned by eigenvectors of , i?2 • • • Rn 
will give reflectivity ~ R' . However, we note that en- 
hanced absorption can also be measured directly, e.g. via 
fluorescence or heat generation. 

We employ a model previously introduced and studied 
by Beenakker et. al. [10] , and by Bruce and Chalker ^ . 
As shown in Fig. [ija), it consists of a waveguide with 
N channels, containing a medium with elastic scattering 
mean free path ^ and ballistic absorption length la- The 
ratio of these lengths defines the absorption parameter 
a = i/^ai with a <C 1 in the regime of interest. In addi- 
tion, the quantity La = gives the average distance a 
diffusing particle penetrates into the medium before be- 
ing absorbed. From the combination law for infinitesimal 
scattering segments (each of which randomly scatters be- 
tween channels and also imposes some absorption param- 
eterized by a), the joint probability distribution (JPD) 
for the reflection eigenvalues can be shown to evolve with 
the sample length L according to the equation [TO] [11] 

'^'^ + ')^ = 2^^ + 2^^i^^^ 

where Px{{xn}) is the JPD given in terms of the variables 
Xn = ^ cosh~^[(l + Rn)/{1 — Rn)]- This change of vari- 
ables leaves the diffusion term in ([T]) independent of {x^}, 
allowing analytic solution of the stationary distribution. 
The "potential" V in Eq. ([T]) is given by 

V = Y^ ^a{N + 1) cosh(2xn) - In | smh{2xn)\ 

n 

— ^ In [cosh(2x^) — cosh(2xn)] . (2) 



Eq. ([T]) is similar to the Dorokhov-Mello-Pereyra- 
Kumar (DMPK) diffusion equation [12^ J^, which is 
derived from a quasi-one-dimensional (quasi-lD) waveg- 
uide model without absorption. The DMPK equation 
describes the JPD of transmission and reflection eigen- 
values, and its results have been found to generalize to 
geometries with no transverse confinement, as when a 
surface is illuminated by a spot [14^ (the geometry used 
in subsequent experiments on coherently enhanced trans- 
mission [7] ) . We expect the results for the absorbing case 
to generalize in the same way. 

For La < L, Px{xir'' i^n) ^ is the stationary 
(L-independent) limiting solution to ([T]). To make con- 
tact with standard random matrix theory (RMT), a fur- 
ther change of variables = 2a{N + l)[cosh(2xn) — 1] 
is performed, which yields 

m>n 

This is the Laguerre eigenvalue distribution, charac- 
teristic of a well-studied ensemble in RMT [15]. The 
ensemble- averaged mean reflectivity can be determined 
by integrating over this distribution JOl [11] : 

{R) ^ l + 2a-2v'a(l + a). (4) 

Note that this result is independent of N. 

We are interested in the distribution of the smallest 
reflection eigenvalue, Ri. This can be found from the 
fact that Nyi follows a xl distribution [15]. Thus, 

-2aiV(Ar+l)^ 

p{R,) = 2aN{N + 1) (i_^^)2 ^ (5) 
and the ensemble- averaged value is 

(Ri) = l + 2a7V(7V + l)e2^^(^+^)Ei[-2aAr(Ar + l)] 
^ [2aN{N + 1)]"^ for a > N'^^^. (6) 

Unlike (R)^ this quantity does depend on N. The de- 
crease of Ri with increasing N is very natural, because 
the eigenvalue repulsion characteristic of RMT ensembles 
tends to push the smallest eigenvalue to small values. 

Thus, even if the material absorptivity is weak enough 
that (R) ~ 1, if A/" is sufficiently large then Ri ~ 1/aN'^ 
can be very small. In typical free-space optical experi- 
ments N ^ A/X^ where A is the spot size, so this condi- 
tion should be easily achieved. According to Q and (|6|, 
the regime of weak average absorption and strong optimal 
absorption is a <C 1 <C 2aN'^. The physical intepretation 
of these inequalities is that photons have little chance of 
being absorbed between scattering events {£a ^ but 
those in the maximally-absorbed mode have a negligible 
chance of diffusing the entire quasi- ID localization length 
^ = Ni before being absorbed. For example, for a < 0.01 
and N = 80, we have {R) > 0.8 and (^i) < 0.01. 
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FIG. 2: (color online) (a) Mean reflectivity R (open sym- 
bols) and smallest reflection eigenvalue Ri (filled symbols) 
vs the absorption strength a. Each data point is an aver- 
age over 2000 disorder realizations. The scattering medium 
has ni = 1.5, disorder strength do = 0.6 (i — 0.05L), and 
uniform absorption. Two simulation sets are shown: (i) ex- 
ternal index no = 1.5 and = 80 exterior waveguide chan- 
nels (circles); and (ii) index mismatch with no = 1 and 
N — 4Q exterior channels (triangles). In each case, we fit 
{R) to Eq. Q to obtain the constant co in — co?7/c/2ni, 
finding co = 2.14 for (i) and co = 4.65 for (ii). This de- 
termines the analytic form of plotted as solid curves, 
with no free parameters. The vertical dashed line denotes 
a > (i/L)^, required for the stationary solution to ap- 
ply, (b) Distributions of Ri and R2 from simulations (his- 
tograms). The solid line shows the analytical expression ^ 
for p{Ri); dotted lines show the ensemble averages (Ri) and 
{R2)' (c) Semi- log plot of reflection eigenvalue density, com- 
paring simulations (histogram) to the approximate density 
2a{N+l)y/[N{l - R)/a{N + 1)R] - l/7r(l-i?)^ (solid line), 
from the large- density of the Laguerre distribution [15]. 
The simulations for (b) and (c) were performed with no = 1.5 
{N = 80) and ry = 0.0003 (4 = 125^). 



To test and generalize the above results, we have per- 
formed numerical simulations of scalar wave propagation 
in quasi-lD structures. The wave equation (V^+e/c^)V^ = 
is discretized using a 2D tight-binding model on a 
square grid. A semi-infinite uniform waveguide of in- 
dex no is connected to a scattering region with average 
index ni, uniform absorption 77 and random scatterers; 
see Fig. [TJa). The reflection matrix is computed for this 
rectangular region, of varying length and width and per- 
fectly reflecting walls on three sides, using the recursive 
Green's function method The elastic mean free 

path i is computed using the relation (T) = i/L^ via sep- 
arate transmission simulations without absorption. In all 
simulations, we use frequency k = 750 /L, lattice spacing 
I// 750, and waveguide width L/5. The ballistic absorp- 
tion length is = cor]k/2ni^ where cq is a constant of 
order unity determined by a fit to Eq. Q. 

The variation of (R) and (Ri) with i/£a is shown in 
Fig. [2|a), together with the analytic predictions. Two 
sets of simulations are plotted. In the first, the lead and 
scattering medium are index-matched (no = ni = 1.5) 



as in the model of [TOl[TT]. In the second, they are mis- 
matched with no = 1, ni = 1.5, representing scattering 
from air into a denser medium, a case for which analytic 
results were not known. In both cases, despite weak aver- 
age absorption {ia > 40^ and (R) ^ 1), we obtain values 
of (Ri) as small as 10~^ to 10~^. For the index-matched 
results the data for (Ri) agree very well with the theory, 
Eq. ([6|, using the number of exterior channels, N = 80. 
For the index-mismatched case, the exterior waveguide 
has fewer channels. A/" = 49, at the same frequency, com- 
pared to the average of 80 interior channels; nonetheless, 
we find that Eq. ([6| holds very well using the exterior 
channel number. Thus, illuminating from a lower index 
medium (e.g. from air) does diminish the CEA effect 
compared to the index-matched case, but large enhance- 
ment is still possible and is governed by the same factor of 
{2aN'^)~^. Figs. [2]^b) and (c) show the numerical distri- 
butions of the reflection eigenvalues for fixed a, which are 
in excellent agreement with the theoretical distributions. 
The distribution of R2 in Fig. |2|b) illustrates the role of 
eigenvalue repulsion in confining Ri to small values. 

The analytic predictions for (Ri) and (R) have no de- 
pendence on the frequency /c, apart from the weak varia- 
tion of the number of exterior channels N with /c, assum- 
ing that the parameter a varies negligibly with frequency. 
Unlike the CPA systems studied in Refs. pL|l2j, in which 
near-perfect absorption is obtained at particular discrete 
values of a and /c, in the present systems highly enhanced 
absorption is possible at any frequency so long as the in- 
cident waveform is optimized in the vicinity of each fre- 
quency of interest. Thus, while CEA is not a broadband 
effect in the usual sense, it can be realized over a wide 
range of frequencies. To determine the frequency interval 
over which a fixed waveform, optimized at a single fre- 
quency, will still lead to strongly enhanced absorption, 
we calculate the frequency correlation of shown in 
Fig. [3j We find that the decorrelation frequency scale is 
kc ^ 0.5/^a7 which is to be expected; by analogy to the 
case of transmission through a non-absorbing medium, 
where the decorrelation scale is determined by the length 
of paths that traverse the sample and escape, here the 
maximally- absorbed mode arises from the interference of 
many coherent paths of length ~ ia. Thus, measure- 
ments obtained at frequency intervals Ak ^ 0.5/ ia are 
independent. For i = 1/im, ia = 100£, and A = 1.5/im, 
this gives AAc ~ 1.5 nm. While the CEA effect is more 
robust than CPA when 1, it should be noted that, 

when A^ ^ 1, continuous CEA is no longer possible, 
whereas CPA can still occur at discrete frequencies. 

An interesting and important extension of the above 
model is to a localized or extended absorber buried be- 
hind a layer of lossless scattering medium. Such a system 
has been studied experimentally by Vellekoop et. a/., who 
demonstrated strong focusing of light through a lossless 
medium to enhance fluorescence in a small interior spot 
[17]. To our knowledge, there has been no theoretical 
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FIG. 3: (color online) Normalized correlation function for 
the smallest reflection eigenvalue defined as CR^{5k) oc 
[{Ri{k + Sk)Ri{k)) - {Ri{k)f] with ^^^(O) = 1. Each data 
set consists of 3000 disorder realizations. Inset: plot showing 
the linear dependence on 1/ia of the decorrelation frequency 
kc, defined by kc = CR^{kc/2) = 1/2. Here, i = 0.05L. 
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interacting gas particles with a stationary distribution at 
length L, in which the confining potential due to a in ([2| 
is balanced by diffusion and interparticle interactions [6j . 
This is the distribution for which {Ri) ^ (2a A/"^)"^ for 
the extremal particle. Adding the lossless region corre- 
sponds to turning off the confining potential for a "time" 
T = Lp/£{N + 1), causing the particles to drift and dif- 
fuse to larger values of Xn- In fact, Fig. Bb) shows that 
{xi{Lp)) - T, so that {Ri) ^ (xj) - T*, which is the 
saturation scale observed in Fig.|4|a). For Lp < i{N-\-l), 
these results show the feasibility of enhancing the delivery 
of energy to a buried absorber by orders of magnitude. 

A final critical question is whether incident waveforms 
approaching the CEA optimum can be found when de- 
tails of the scattering medium are unknown, via an em- 
pirical optimization scheme analogous to those employed 
for transmission through lossless media ^ (8^ . Initial nu- 
merical tests based on our model are promising, with 
simple optimization schemes giving order of magnitude 
increases in absorption, though not yet approaching the- 
oretical limits. This problem, and that of the buried ab- 
sorber, are rich and exciting directions for future work. 
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FIG. 4: (color online) (a) Variation of Ri with a for buried 
absorbers. In the first set of simulations (triangles), the ab- 
sorbing region (length L, width 0.2L equal to the waveguide 
width, and i = 0.05L) lies behind a segment of lossless scatter- 
ing material (length Lp — lOi, with the same i). In the second 
set of simulations (crosses), the absorbing region has dimen- 
sions 0.06L X 0.08L and is surrounded by the lossless medium. 
The uniform absorption case, similar to Fig. [2] is shown for 
comparison (circles), (b) Plot of ^y (Ri) vs the Fokker-Planck 
"time" Lp/£{N-\-l), for the full- width absorber with a = 0.04. 
Note that Rn ~ x^, where Xfi are the variables of Eq. M. 

work on the limiting efficiency of such a process. We con- 
sider two variant models, shown in the inset of Fig.|4ja). 
A lossless region of length Lp is added in front of the 
absorbing region of length L, which contains either uni- 
form absorption or a localized absorber. For a within 
the absorber satisfying 2aN'^ ^ 1, the simulations still 
show GEA of approximately two orders of magnitude. 
However, Ri now saturates with increasing a to a value 
« (Lp/iN)'^. This saturation value appears to be inde- 
pendent of the size of the absorbing region, though the 
larger region reaches this value for smaller a. 

To understand this behavior qualitatively, we return 
to Eq. ([T])-([2|. One can think of the eigenvalues {xn} as 
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